Predictions of the causal entropic principle 
for environmental conditions of the universe 
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The causal entropic principle has been proposed as a superior alternative to the anthropic principle 
for understanding the magnitude of the cosmological constant. In this approach, the probability to 
create observers is assumed to be proportional to the entropy production AS in a maximal causally 
connected region — the causal diamond. We improve on the original treatment by better quantifying 
the entropy production due to stars, using an analytic model for the star formation history which 
accurately accounts for changes in cosmological parameters. We calculate the dependence of AS on 
the density contrast Q = Sp/p, and find that our universe is much closer to the most probable value 
of Q than in the usual anthropic approach and that probabilities are relatively weakly dependent 
on this amplitude. In addition, we make first estimates of the dependence of AS on the baryon 
fraction and overall matter abundance. Finally, we also explore the possibility that decays of dark 
matter, suggested by various observed gamma ray excesses, might produce a comparable amount of 
entropy to stars. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

The nature of the cosmological constant problem has 
changed since it was fully recognized in the 1980s. At 
that time, a seemingly natural solution was to find a 
mechanism which would explain why A should be zero. 
Among the most promising such ideas were those of 
Hawking [l[ , elaborated by Coleman Q , in which the Eu- 
clidean path integral for gravity was interpreted as a wave 
function for the universe, shown to be strongly peaked at 
A = 0. However we now know with a high degree of con- 
fidence that A is not zero, due to observations of distant 
type la supernovae combined with constraints on the 
flatness of the universe from the CMB [i| and comple- 
mentary information from the x-ray baryon fraction in 
galactic clusters [5| . 

Two notable ideas from the 1980's predicted that A 
should be nonzero: the membrane creation mechanism 
of Brown and Teitelboim (BT) [|| (building on work in 
0), involving tunneling between false vacua with differ- 
ent values of A, and the anthropic argument of Weinberg 
0, which predicted that 



A obs < A<100A obs 



(1) 



based on the requirement that galaxies should be able to 
form and thus give rise to observers before the universe 
collapsed or expanded too quickly. These approaches 
have recently gained support in a number of ways from 
the string theory community: by the explicit realization 
of the BT mechanism within heterotic M-theory and type 
II string theory [9f and the realization that the vacuum 



'Electronic address: jcline@physics. mcgill.ca 
^Electronic address: frey@hep. physics. mcgill.ca 
* Electronic address: holderOhep. physics. mcgill.ca 



structure of string theory is a vast landscape llOl e ncom- 
passmg more than 10 500 possible values of A [ll]. The 
picture that has emerged is that this multitude of vacua 
can be simultaneously realized through the process of 
eternal inflation (T^ |. Tunneling a la BT from the eter- 
nally inflating regions would populate universes like ours 
with small A, as well as those with much larger A. The 
latter would be incompatible with life, and we should dis- 
regard them if we are interested in explaining properties 
of universes which can admit the existence of observers. 
This is the weakest form of the anthropic principle; it 
does not insist that physics must be compatible with our 
existence, only that we take our existence as a data point 
in putting experimental constraints on the landscape of 
vacua. 

This point of view has been strongly criticized on the 
grounds that it admits no real predictions, only post- 
dictions. Moreover, the anthropic explanation of A has 
been weakened by the discovery of dwarf galaxies at 
z ~ 10, proving that structure formed earlier than z ~ 4 
as was supposed in [8j. Since the bound (JTJ scales like 
(1 + z) 3 , the anthropic constraint is weakened by an or- 
der of magnitude. To counteract this, it has been argued 
that only galaxies of a certain minimum size are capa- 
ble of retaining the heavy elements needed to sustain life 
[l3| . but such assumptions seem questionable and make 
the anthropic principle appear to be increasingly arbi- 
trary. The anthropic bound on A is further weakened 
by allowing other cosmological parameters to vary at the 
same time. For example, the primordial density contrast 
Q = Sp/p ~ 2 x 10 -5 in our universe, but, if Q were 
larger, structure could form earlier despite a larger ex- 
pansion rate. The anthropic bound on A scales like Q 3 , 
so increasin g Q can relax the bound by many orders of 
magnitude |14T [HI ]. Similarly, [lj[ showed that allow- 
ing the effective Planck scale of a tensor-scalar theory of 
gravity to vary also weakens the anthropic bound on A; 
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this result is potentially important for the string theory 
landscape. 

Despite these shortcomings, one should keep in mind 
that the anthropic approach was the only one to predict 
the range of A before its nonzero observation, and the 
value A obs ~ 10~ 123 M 2 is so peculiar from the parti- 
cle physics perspective that it seems exceedingly unlikely 
that a dynamical mechanism could by itself explain the 
observed value of A. Something like the anthropic prin- 
ciple therefore appears to be necessary for understanding 
the magnitude of the dark energy. One is thus motivated 
to search for some improvement which retains the virtues 
of the anthropic principle while getting rid of overly spe- 
cific assumptions about the nature of observers. Such 
an idea, dubbed the "causal entropic principle," has re- 
cently been proposed by [l7|. It assumes that a good 
tracer of the potential for forming structure, and hence 
observers, is the amount of entropy produced within a 
causally connected region (the causal diamond) defined 
to start at some early initial time £j, such as reheating, 
and ending in the infinite future, which corresponds to a 
finite conformal time r if A > 0. The created entropy, 
AS, does not include entropy already present at the mo- 
ment of reheating but only that which is created after 

k. 

The causal entropic principle has several virtues which 
make it worthy of further investigation. First, it predicts 
a probability distribution dP/dlogA such that A Q b s is 
within 1 a of the most likely value of A, unlike the an- 
thropic approach for which P(A < A Q b s ) < 10 -3 . Sec- 
ond, it makes minimal assumptions about the detailed 
nature of observers or life. It only assumes that free en- 
ergy (leading to increased entropy) is available, which is 
a fundamental thcrmodynamical requirement for making 
any measurement. Third, it provides a specific proposal 
to separate potentially divergent volume factors (for ex- 
ample, from eternal inflation) into the prior probability 
distribution. (For recent discussions of volume-weighted 
measures, see fl8j.) 

In their paper, the authors of [17| understandably re- 
stricted their attention to the determination of A. To 
carry the idea further, we are interested in exploring the 
probability distributions for other cosmological param- 
eters, for example Q. This requires knowing how the 
rate of entropy production depends on these parameters. 
In particular, since [l7| found that starlight is the main 
source of entropy, a major result of the present work is to 
show how the rate of star formation depends on various 
cosmological parameters. This was somewhat roughly 
estimated in [17| : here we have tried to be more quan- 
titative. The most straightforward quantity to vary in 
fact is Q, and we have studied this dependence to show 
that the predicted value of A in the causal entropic ap- 
proach is much less sensitive to Q than in the anthropic 
approach. Furthermore, the entropic prediction for Q is 
much closer to the observed value than is the anthropic 
prediction. Similarly, we find that the prediction for the 
baryon fraction is within a factor of a few of the observed 



value. 

An important aspect of [TtJ is the realization that the 
typical timescale of star formation and galaxy evolution is 
of the same order of magnitude as the age of the universe, 
providing another coincidence problem. The causal en- 
tropic view effectively elevates this coincidence to a prin- 
ciple, with potentially testable predictions. 1 A process 
that generates a large amount of entropy on a character- 
istic timescale would make more likely, by this principle, 
a cosmological constant that becomes dynamically im- 
portant on a comparable timescale. This suggests that 
there is not a large source of entropy that is acting on 
a timescale significantly different than the current age 
of the universe. If we discover a new important entropy 
source, we should expect its timescale to be the cosmo- 
logical one. In fact, such an entropy source would make 
the observed value of A more likely. 

Beyond the consideration of universes very similar to 
our own, the entropic approach has the capability of 
quantifying the likelihood of universes with quite differ- 
ent properties. This leads one to ask whether alternative 
universes with much greater entropy production than our 
own can be constructed, addressing the critique of [2l[ 
that existing work has focused too much on universes and 
life similar to our own. 2 If so, it might call into ques- 
tion the validity of entropic principle, since our universe 
would then appear to be relatively unlikely One way in 
which entropy might be copiously generated is through 
the decay of dark matter. Interestingly, hints of excess 
gamma rays at various energies has led to speculation 
that the dark matter in our universe indeed is unstable. 
One goal of this work is to determine whether the en- 
tropy released by such decays can be competitive with 
that due to starlight; if so, this could be an indication 
that the entropic principle can explain properties of dark 
matter as well as dark energy. 

We will review the causal entropic approach in section 
HI] including our results for dP/d log A, and motivate the 
improved star formation rate (SFR) which we develop 
in section Mil This will be followed in IIVI by our re- 
sults for the probability distribution for Q and for the 
baryon fraction. Section [V] will discuss the production 
of entropy from decaying dark matter. We give conclu- 
sions in section fVTl Throughout, except when units are 
given explicitly, we work with h = c = 1. In formulae, 
we keep factors of Mp and G explicit, but we omit them 



1 An interesting paper [lfj has recently analyzed the joint prob- 
ability of observing a universe with a given A at a given CMB 
temperature using star formation and evolution as the key in- 
dicator of observers, following the original ideas of 12(1. Even 
though the spirit of |19( is not as general as that of |17|. the 
results are, unsurprisingly, similar. 

2 Also, like [2l| . one might ask whether we should be interested in 
the conditional probability of observing A b 8 given the existence 
of observers or given the existence of human-like life. While 
more traditional approaches focus on human-like life, the causal 
entropic principle allows us to consider more generic observers. 
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for brevity in figures and when writing dP/d log A (which 
should be read as dP/dlog(A/M],)). 



II. THE CAUSAL ENTROPIC PRINCIPLE 

In this section, we briefly review the basic ideas and 
results of [T^. The fundamental assumption is that the 
probability distribution for A is proportional to entropy 
production in the causal diamond. 



dP 
dA 



cx AS 



dtV c {t) 



dS 
dV r dt 



(2) 



where V c = (47r/3)r(i) 3 is the comoving volume in the 
diamond at a given time, and dS/dV c dt is the rate of en- 
tropy production per comoving volume. The causal dia- 
mond is defined in conformal time, ds 2 — a 2 (T)(—dT 2 + 
dx 2 ), as the region shown in figure [TJ The initial time 
is taken to be the time of reheating, although results are 
quite insensitive to this choice since V c is negligible at 
such early times. In conformal time, the final time is fi- 
nite (in fact, it is defined to be 0) since r = J dt/a(t) 
converges for any A > 0; recall that a(t) ~ exp(i-^/A/3) 

1/3 

asymptotically. The initial time is — r max rs — 2.8 i A = 
2.8 (3/A) 1 / 6 . The diamond achieves a maximum volume 
at the intermediate time — r max . 




max 



x 



FIG. 1: The causal diamond 



Ref. [17| estimated the amount of entropy produced 
by various astrophysical sources, including active galac- 
tic nuclei, supernovae, and cooling of galaxies, and de- 
termined that starlight, inelastically scattered by dust, 
is the largest source of entropy. 3 The rate of entropy 
production per comoving volume can be written as 



dS 
dVrdt 



(t) = f df 
Jo 



d 2 S 
dM+dt 



(t-t')Mt') , 



(3) 



3 The horizon entropy of black holes is a special case which is 
excluded; a single black hole's entropy far outnumbers that of 
other entropy sources in the universe, but a universe filled with 
black holes would seem to be inimical to observers. 



where (d 2 S/dM+dt)(t — t') is the rate of entropy produc- 
tion per stellar mass at time t due to stars born at time 
t', while p*(t') is the rate of stellar mass production (star 
formation rate, or SFR) at time t'. (d 2 S/dM*dt)(t - t') 
in turn is given by 



d 2 S 
dM+dt 



(t-t') = 



1 



(M) 



<(*-*') 



O.O8M 



dM 6mf(M) 



d 2 t 



dN„dt 



(4) 

in terms of the initial mass function £imf(-^0 (equal 
to 0.105 M~ 2 - 35 for M > 0.5M Q and 0.189 M" 1 - 5 for 
M < 0.5M Q ). The entropy production rate for a sin- 
gle star, d 2 s/dN±dt, is its luminosity over an effective 
temperature, L±/T cS w (M/M ) 3 5 x 10 54 yr" 1 . The ef- 
fective temperature T e g w 20 meV reflects the fact that 
half of all starlight, originally UV, is scattered into the 
IR by dust and its associated entropy is thereby mul- 
tiplied. The integration limits in ^ are the minimum 
and maximum masses of stars of a given age, where 
Mmax(i-t') = max(lOOM , [10 Gyr/(i-i')] 2/5 ), reflect- 
ing the fact that the largest stars burn out more quickly 
than the smallest ones. 

The total entropy production is thus the convolution 
of a cosmological factor, V c (t), with an astrophysical one, 
(dS/dV c dt)(t). The former gives a preference to small 
values of A since V c (t) cx 1/VA and /~ dtV c (t) cx 1/A. 
The latter function is peaked at a time 2 — 3 Gyr, de- 
pending on the choice of the SFR, p*, of which there are 
several somewhat different determinations in the litera- 
ture [13, HI) HI]]- The fact that star formation tapers 
off in the late universe softens the preference for small A 
such that the total entropy production (and hence by as- 
sumption the probability distribution for A) approaches 
a constant as A — > 0, but falls off rapidly as A becomes 
much greater than the observed value. To better visual- 
ize the range of likely A values, it is useful to consider 
the distribution in log A, dP/d\ogA = AdP/dA, which 
is peaked at a nonvanishing value of A. We have repro- 
duced (and extended) the calculations of [ljj to obtain 
dP/dlogA, using three different determinations for the 
SFR, Nagamine et al. j|22] (N), Hopkins-Beacom HU 
(HB), and Hernquist-Springel (HS) [24[. The distribu- 
tions (not normalized) are shown in fig. [5] 

Although there is some uncertainty in the prediction 
for dPj d log A due to differences in the estimates for the 
SFR, all three curves in fig. [5] show that the observed 
value of A is within 1 a of the most likely value, which 
is a much better agreement than that obtained using the 
anthropic approach. In [llj , only the the HB and N SFRs 
were considered, which are phenomenologically derived 
for our universe and in particular for the observed value of 
A. In order to estimate the effect of varying A, [17| multi- 
plied the entire SFR by a correction factor F(t, M), taken 
to be the Press-Schechter fraction 25j of matter collapsed 
into halos at a typical galactic mass scale 10 7 M Q and the 
time t ~ 2 — 3 Gyr when the SFR is maximized. (We 
discuss the Press-Schechter formalism in some detail in 
the appendix.) This is a reasonable, though somewhat 
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FIG. 2: Unnormalized probability distributions for A using 
different star formation rates 22. [23l. [24T|. Vertical line indi- 
cates the observed value of A. 



crude approximation, since the true answer must depend 
on the details of gravitational collapse over a range of 
mass and time scales. A major goal of the present work 
is to improve on this approximation by developing an 
SFR which has the correct quantitative dependences on 
A and other cosmological parameters. We have carried 
this out, based on an analytical model due to Hernquist 
and Springel [H| (henceforth HS), which is the subject of 
the next section. Our result for dP/d log A based on the 
HS SFR, as shown in fig. [2j is in reasonable agreement 
with the earlier results. More importantly, it allows us 
to consider with greater confidence the probability distri- 
butions for other parameters, as a way of further testing 
the causal entropic principle. 



III. THE HERNQUIST-SPRINGEL SFR 

As mentioned above, the N [HJ and HB [H[ SFRs 
are simple phenomenological formulae derived from ob- 
servations in our universe, so they cannot tell us how the 
variation of cosmological parameters will affect the SFR. 
Fortunately, Hernquist & Springel [lij have developed a 
simple analytic model for the SFR based on more de- 
tailed numerical simulations [26[ . Because the HS model 
can be written explicitly in terms of cosmological param- 
eters and other physical constants, it is possible to vary 
the SFR in reponse to changes in, for example, the cos- 
mological constant or the amplitude of cosmological per- 
turbations. In this section, we will briefly describe our 
implementation of the HS SFR model, leaving a more 
detailed discussion of the model for the appendix. 

HS start with the formula 



dF 
rflnM 



(M,t) s(M,t) rflnM . 



(5) 



of matter collapsed into clouds of mass M or less) [25| , 4 
and s = (f>i,) is the averaged rate of star formation in 
collapsed haloes. 

The collapsed fraction for a given scale as a function of 
time can be simply calculated from the statistics of Gaus- 
sian random fields for an assumed input matter power 
spectrum with a power spectrum amplitude that evolves 
with time as expected from linear theory. 

The scale at which the collapsed fraction should be 
calculated is set by atomic physics: one can collisionally 
excite hydrogen atomic transitions when kinetic temper- 
atures exceed ~ 10 4 K. The virial relation between mass 
and temperature at time t is found to be 



M = 



(2fcT//i m ) 3 / 2 _ (2fcr/ Mm ) 3 / 2 
p v ~ lQGH{t) 



(6) 



where k is Boltzmann's constant, fi m is the mean mass 
per particle, and M depends on the Hubble parameter 
at the given time. 

Using constants S c — 1.6868 and a — 0.707, num- 
bers found to provide good agreement beween theory and 
large N-body simulations for the statistics of collapse, a 
simple version of the SFR suggested by HS is 



(7) 



where 174 stands for fluctuations at the mass scale which 
virializes at temperature 10 K, and the physics of cooling 
is parameterized in the prefactor q(t). 

The star formation rate within a collapse object is reg- 
ulated by the rate of radiative cooling and the efficiency 
of radiative cooling is difficult to calculate; HS report 
that a good approximation to the star formation effi- 
ciency can be characterized by 



q(t)^[ X (t)x/(x(t) m + X m ) 1/m } 9/2V 



(8) 



where x = (H/Hq) 2 / 3 , and x = 4-6, 77 = 1.65, and m = 6 
provide good fits to numerical simulations (Ho is the 
observed value of the Hubble rate taken as a reference 
value). Because we are interested in relative probabili- 
ties, we do not need to normalize the star formation rate. 

To illustrate the parametric dependences of our mod- 
ified SFR, we plot p*(i) for several values of A in figure 
131 We also show there a naive approximation to the HS 
result, which consists of the analytic fit which HS made 
to their result, 



c(t) <x 



X 



(l + a( X -l) 3 exp(/V/4) 



(9) 



Here, F(M, t) is the Press-Shechter fraction (the fraction 



4 HS advocate the use of a more precise collapse fraction derived by 
Sheth & Torman [27|, but the Press-Schechter fraction is accurate 
enough for our purposes. 
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Here all dependence on the cosmological parameters is 
due to x, which depends upon A through the time- 
dependent Hubble rate, and the parameter /3, which de- 
pends on the density contrast Q as j3 ~ Q~ 2 via the 
Press-Schechter formalism. We have normalized © to 
agree with the more exact result at late times. The rough 
approximation tends to overestimate or underestimate 
the peak value of the more accurate SFR. 

Let us also contrast our treatment with that of [l~7l |. 
where the shape of the SFR was assumed to remain con- 
stant, and only its overall normalization depended on 
A. Our results indicate a significant change in shape 
with A, with the peak even disappearing at large A. The 
reader may be surprised that the SFR increases with A, 
since increasing the expansion rate is supposed to delay 
formation of structure. However since the effect of A 
appears in the prefactor of p*, through x ( see ec l s - ® 
or ©), there is initially an increase with A. The de- 
crease predicted by Press-Schechter becomes apparent 
for log(A/A/|,) > — 121. In figure H] we show this de- 
pendence by plotting the large-time asymptotic value of 
p\ as a function of A. This function is peaked around 
log(A/Mj.) > - 121. Taken by itself it would indicate 
a lower probability for our universe with log(A/Af|>) ~ 
— 123 than the full entropic treatment gives. 




time (G yr) 

FIG. 3: Solid lines: the star formation rate used in this paper, 
for several values of p\- Dashed lines are the approximation 
(J9]). Common proportionality constants have been omitted. 

In the above we have given the simplest reasonable 
implementation of the HS model, and there is room for 
improvement, as outlined in the appendix. However, this 
is a physically motivated model for how star formation 
depends on cosmological parameters, and it captures the 
important elements that any such model must have. 

IV. PREDICTIONS FOR OTHER 
PARAMETERS 

As long as the properties of the universe are such 
that starlight continues to be the dominant source of 
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FIG. 4: Asymptotic value of the star formation rate at large 
times, as a function of A/M|>. 



entropy production, we now have the tools to compute 
AS as a function of cosmological parameters beyond just 
A. These include the density contrast Q, the ratio of 
baryonic matter to dark matter, the ratio of matter to 
photons, and possibly the spatial curvature. The lat- 
ter is more subtle to try to quantify because it is time- 
dependent. One must specify the curvature at some 
reference time, for example at matter-radiation equal- 
ity. However, there are good reasons for believing that 
the curvature of our universe was determined by infla- 
tion rather than by environmental considerations. We 
therefore will confine our investigations to the explicitly 
time-independent quantities. 



A. Density contrast 

The most straightforward cosmological parameter to 
vary is the density contrast Q, since it appears only 
through (|A.24|) . It is also interesting because of the 
great sensitivity of the anthropic prediction for A on Q. 
This can be understood from (|A.10p and (|A.24|1 . which 
show that Q and A appear in the Press-Schechter frac- 
tion through the combination QA 1//3 . In the anthropic 
approach, structure formation is the most important ef- 
fect, and this causes the anthropic upper bound on A 
to scale like Q 3 . In contrast, the entropic approach gives 
less weight to very early structure formation due to larger 
Q, because the volume of the causal diamond is smaller 
at earlier times. One can thus anticipate less sensitivity 
to Q in the entropic approach. 

This expectation is borne out in figure [5J where Q is 
varied by 4 orders of magnitude above the observed value, 
Q = QsQobs- In fact, the curves for dP/dlogA saturate 
near the highest one shown as Q — + oo. The position 
of the peak of the distribution shifts by only an order of 
magnitude, which is a much weaker dependence than the 
A ~ Q 3 dependence due to structure formation alone. 
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FIG. 5: Unnormalized probability distributions for A at dif- 
ferent values of Q = Q S Q obs 
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FIG. 6: Unnormalized probability distribution for Q at fixed 
A = A obs . 



As Q is varied, the epoch of star formation moves to 
earlier times, but the characteristic time of 1 Gyr remains 
relevant. It matters little whether stars form at 10 6 years 
or 10 8 years since their characteristic lifetime is 10 9 years. 
At extremely low Q, though, the A ~ Q 3 dependence 
should reemerge. 

We also plot AS(Q) for fixed A = A b s in figure [6] and a 
contour plot varying both A and Q in figure [71 Although 
our universe is not at the peak of the distribution, the 
probability of Q b s is only a factor of 8 smaller than the 
most likely value. For comparison, the usual x 2 statis- 
tic is of order — 21nP, so a factor of 8 corresponds to a 
change in the effective \ 2 on the order of 4.2. 

One effect which has not been quantified in our ap- 
proach is the cutoff on Q which would result from over- 
production of black holes; beyond some large value of Q, 
formation of stars (and observers) would be impeded by 
the loss of material to black holes. A discussion of this 
issue can be found in 

The reader should note that we have only calculated 
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FIG. 7: Contours of the probability distribution dP/d log A, 
varying Q and A. 



the entropic weighting factor for the probability, inde- 
pendent of any prior distribution for Q. What our re- 
sults show is that, first, the observed value of A is still 
near the most likely value, even for large values of Q, 
and also that entropic effects will not modify the prior 
distribution for Q by more than an order of magnitude 
across a large range. It is reasonable to expect the prior 
distribution for Q to be dP/dQ cx Q n with —l^n<l 
[191; the smaller value works well with our results. 



B. Baryon fraction 

Another cosmological parameter which can obviously 
affect star formation is the ratio of baryonic to dark mat- 
ter, or baryon fraction fa. Since stars are baryonic, the 
SFR should have an overall factor of fa to account for 
the availability of baryonic material for star formation, 
as appears in (|A.16|1 . 

The details of structure formation are also affected by 
fb] at very large values, gravitational collapse is impeded 
by pressure; this is known as Silk damping. Because the 
baryonic matter cannot collapse easily due to pressure, 
the matter power spectrum in (|A.4j) scales as (1 — /&) for 
a fixed initial amplitude Q, so the r.m.s. fluctuation a 
scales like (1 - fa) 1 / 2 . 

Since our universe has fa ~ 1/6, wc define the relative 
baryon fraction as r& = 6 fa. The dependence of on 
rft as described above can be summarized as making the 
replacement 



p*(jj,, cr) -» r b p* /I, 



6 - n 



1/2 



(10) 



The entropic distribution for is shown in figure[5J Sim- 
ilarly to the distributions for A and Q, it indicates that 
the measured value of fa is not quite the optimal one, 
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but also not very far from being so: the probability of 
/b.obs is ~ 1/2 of that for the most likely value. 

The reader should note that, as with Q, we have not 
attempted to calculate the prior probability distribution 
for the baryon fraction. Our results show that the en- 
tropic weighting does not change the probability distri- 
bution by more than a factor of 2 for 1/6 < /& < 5/6. 
Therefore, the value taken in our universe is reasonably 
probable unless the prior distribution is strongly peaked 
at some other value. 

25 r 1 i 1 i 1 i 1 i 1 i 1 =i 
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relative baryon fraction 



FIG. 8: Unnormalized probability distribution for relative 
baryon fraction r b = 6/f, at fixed A = A obs . 



£ increases but in a manner qualitatively similar to the 
dependence on the density contrast Q. This is shown in 
figure [HI Similarly to the case of Q, the relative proba- 
bility of the observed value of £ is ~ 1/8. 

Again, we have not attempted to calculate the prior 
probability distribution for £, but our results show that 
the entropic weighting factor itself does not disfavor our 
universe greatly. 
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FIG. 9: Unnormalized probability distribution for relative 
matter abundance r m at fixed A = A bs- 



C. Matter abundance 

A third quantity which can strongly affect star for- 
mation is the overall abundance £ of matter relative to 
photons, for fixed baryon fraction, Q and A. In a uni- 
verse with larger £, matter domination will occur ear- 
lier, and the CMB temperature at a given time will be 
reduced. In our analysis, radiation does not play a sig- 
nificant role in the evolution of the scale factor of the 
universe, so we can again focus on how £ affects the SFR. 
The dependence of the SFR on £ has two origins; first, 
fi scales as £ 2 by definition. Second, the linear growth 
factor for perturbations, G(a eq ) in (|A.8j) . depends on £ 
simply because it is evaluated at matter-radiation equal- 
ity. As long as matter dominates over curvature and A 
at matter-radiation equality, it is straightforward to find 
[lU that G(a eq ) ~ £ -4 / 3 . From those arguments, we can 
infer that the £ dependence comes in the form 

ct ) P* ( r ™^' r rl 3(T ) > ( n ) 

where the relative matter abundance r m is the rescal- 
ing of £ such that r m = 1 for our universe. Unlike the 
case of varying /;,, there is no reduction in the matter 
power spectrum as r m is increased, so star formation 
and entropy production only become more efficient as 



V. DECAYING DARK MATTER 

In the previous section we confined our attention to 
the variation of standard cosmological parameters. More 
generally, one could imagine universes with properties 
rather different from our own, requiring for their descrip- 
tion other parameters than those characterizing the stan- 
dard cosmological model. A possible new source of en- 
tropy is the decay of dark matter particles. We shall 
show in this section that one can design universes not 
so different from our own where entropy production is 
actually dominated by dark matter decays. 

Interestingly, there are experimental hints from obser- 
vations of gamma rays that dark matter could be unsta- 
ble. The EGRET collaboration observes 7-rays in the 
2 — 10 GeV range which are unaccounted for by standard 
mechanisms [33| . It has been suggested that these excess 
7s are due to the annihilation [34[ or decays [35[ of dark 
matter. At lower energies, the COMPTEL experiment 
observed an excess in the 1 — 5 MeV range. Ref. [36| 
has proposed the decays of Kaluza-Klein (KK) or super- 
symmetric (SUSY) dark matter to explain this anomaly. 
In a similar energy range, numerous experiments, includ- 
ing the SPI spectrometer on the INTEGRAL observatory 
37j have observed 511 keV photons from the galactic 
center, indicative of an excess of positrons annihilating 
nearly at rest. Astrophysical explanations have not con- 
vincingly accounted for this excess, leading to suggestions 
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that the positrons are a result of dark matter annihila- 
tions or decays (39| . 

If indeed dark matter is shown to be unstable, it is 
curious that the MeV-scale anomaly requires a lifetime 
which is within a few orders of magnitude of the age of 
the universe. This could then pose a new coincidence 
problem. Interestingly, the entropic principle could ex- 
plain such a coincidence, as we now discuss. The rate 
of entropy production from dark matter (DM) decays is 
much simpler to calculate than that from stars; it is given 
by 



dS d _ rt 
= g s l ne 



dVrdt 



(12) 



where T is the DM decay rate, n the DM density, and 
g s the entropy increase per particle decay. We find an- 
alytically that the entropy produced by decays has the 
functional form 



AS d = t A f(Tt A ) 



(13) 



where we recall that t A = y/3/A. The function /, calcu- 
lated numerically, is shown in figure 1101 It is closely fit 
by the analytic approximation 



In / w ci + c 2 x - ln(l + c 3 e c 



(14) 



where x = lnl^A and the constants are c\ — 1.164, C2 = 
0.0512, c 3 = 0.889, and c 4 = 1.921. The distribution is 
peaked for a DM lifetime given by 



= I ~ ±L 

T r ~ 4.7 



(15) 



In our universe, t A = 16.7 Gyr, corresponding to a pre- 
ferred lifetime of 3.6 Gyr. This could then explain why 
dark matter should be decaying on a timescale compara- 
ble to the age of the universe. 




FIG. 10: Entropy produced by dark matter decays as a func- 
tion of the dimensionless combination TtA (upto an overall 
proportionality) . 



coming from starlight. We must therefore compare the 
magnitudes of AS from decays and from stars. Ref. [ItJ 
found that the peak rate of entropy production by stars 
was of the order 

' /S * lQ^Mpc-'y" 1 ■ (16) 



dVrdt 



Using the known energy density of dark matter, we find 
that the corresponding expression from decays is 



dS d 
dVrdt 



10 66 3 . 



( r 



VGyr" 



Pdm 



eV\ 

^Jv( 10 ~ 3eV )' 



Mpc _; y 



where m is the DM mass. Therefore the requirement for 
entropy from DM to dominate is 



171 , , r 

— < keV 

9s 



(18) 



This can obviously be satisfied even if g s ~ 1 if dark 
matter is sufficiently light. 5 However, to explain the 511 
keV signal, for example, one needs m ~ MeV, which 
would require g s ~ 1000. At first sight, this would seem 
like an unreasonably large increase in entropy per de- 
cay. However, it might actually be easy to achieve, when 
one considers that electrons are typically produced along 
with the positrons, since we know that dark matter is 
charge-neutral, and that these electrons, if even mildly 
relativistic, must produce a great number of lower-energy 
photons as they thermalize. 

Let us consider how an electron of energy ~ MeV ther- 
malizes within a galaxy. Synchrotron radiation is kine- 
matically blocked because the cyclotron frequency is only 
eB/m e ~ 10" 13 eV for B ~ 10~ 9 T, while photons have 
a plasma mass of m 7 ~ (ane/nie) 1 ^ 2 ~ 10~ 12 eV, since 
n e ~ 1 cm -3 . Instead, the electron loses energy to the 
galactic medium by Coulomb interactions; using the for- 
malism developed by we can see that an electron 
would lose all its kinetic energy over ~ 10 5 yr and ~ 3 
kpc. Thermalization thus takes place much faster than 
the age of the universe. The final energy is of the or- 
der E r*> 1 — 100 eV, resulting in an entropy production 
of AS ~ MeV/E ~ 10 4 - 10 6 . For electrons produced 
in the galactic halo, thermalization will be less efficient. 
Nevertheless, from this number we see that it is possible 
for DM decays which produce mildly relativistic e + e~ 
pairs to result in entropy production which is similar in 
magnitude to that produced by stars. 

It might be objected that the existence of such low en- 
ergy photons should be irrelevant for creating observers, 
even if they far outnumber photons originating from 
stars. However we wish to avoid making any assump- 
tions about the detailed nature of observers, since this 
was one of the arbitrary features of the anthropic princi- 
ple which one would like to overcome with the entropic 
principle. 



The above argument might only be valid if the en- 
tropy produced by decays is greater than or equal to that 



J Notice that g s is quantized, so it cannot be arbitrarily small and 
still nonzero. 
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VI. CONCLUSIONS 

Using the analytic star formation model of HS, we have 
improved upon and extended the calculations of to 
test the causal entropic principle further. First, we recal- 
culated the probability distribution for the cosmological 
constant A alone, finding good agreement with the ap- 
proximations of [17( . We subsequently extended the orig- 
inal analysis to allow for simultaneous variation of A with 
the density perturbation amplitude Q, finding that the 
most probable value of A changes by less than an order 
of magnitude as Q varies by 5 orders of magnitude. In 
addition, the width of the probability distribution grows 
as Q grows, so the observed value of A remains reason- 
ably probable even at large Q. The observed value of Q 
(at fixed A) has about 1/8 the probability of the peak 
values, which may or may not be compensated by the 
prior probability distribution. Due to the flexibility of 
the HS SFR, we have also been able to calculate the en- 
tropic weighting factor for varying baryon fraction and 
matter abundance, finding that the observed values are 
not unreasonable. 

We have moreoever demonstrated the adaptability of 
the causal entropic principle to universes which could be 
qualitatively different from our own, by considering en- 
tropy produced by the decay of massive particles. If en- 
tropy due to particle decay is the dominant source of en- 
tropy production, we find that the particle lifetime should 
be about 20% of t\ (i.e. 3.6 Gyr, since t& = 16.7 Gyr 
in our universe) for the optimal increase in entropy. We 
showed that in fact it is easy for entropy from decays to 
dominate that from stars if the dark matter is warm (m ~ 
keV) . Although warm dark matter seems to be ruled out 
in our universe, a small warm component could still gen- 
erate significant entropy through its decays. Intriguingly, 
even decaying MeV dark matter, hinted at by various 
observed excesses in the galactic and diffuse gamma ray 
spectra, could produce more entropy than stars if its life- 
time is comparable to the age of the universe. 

In summary, we have tested the causal entropic prin- 
ciple by varying a number of parameters which could be 
environmental in nature, as opposed to being fundamen- 
tal constants, and it has so far survived these tests. More 
ambitiously, one could consider particle physics param- 
eters like the masses and charges of electrons and pro- 
tons varying in a landscape of vacua of the fundamental 
underlying theory. Even further afield, we can imagine 
universes with different gauge groups and particle con- 
tent. It may be interesting to consider how the entropy 
produced in such universes compares to our own, to get 
a better idea of whether the entropic principle works in 
a broader setting. 
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APPENDIX: DETAILS OF THE HERNQUIST & 
SPRINGEL SFR 

In this appendix, we review in more detail the HS SFR, 
which was developed based on more detailed numerical 
simulations [26j combined with analytic reasoning. 

We remind the reader of the basic formula 

f dF 

P*(t) = / j^(M,t)s(M,t)dlnM . (A.l) 

Here, F(M, t) is the Press-Shechter fraction (the fraction 
of matter collapsed into clouds of mass M or less) [25[ , 
and s = (p*) is the averaged rate of star formation in 
collapsed haloes. We can now break down the individual 
parts of this formula. 



1. Virialization 

An important issue for star formation is the virializa- 
tion of collapsing gas, which we here review, following 
the discussion of [28j on the spherical collapse model in 
a matter-dominated FRW universe. 

Consider a universe of critical density in matter with 
a spherical overdensity (and a corresponding underdense 
shell). There is an exact solution for the evolution of 
this universe starting from a = at proper time t = 
(its precise form is unimportant for us). The overdense 
sphere reaches its turn-around radius at time t ta and then 
recollapses completely (to a local value a = 0) in time 
2t ta - At the time of complete recollapse, the overdensity 
in the linearized theory is S c = 1.686, which is known 
as the critical overdensity. Any fluctuation in the initial 
density field with a larger linearly evolved overdensity 
would have collapsed at an earlier time. 

Of course, the overdense sphere does not recollapse 
completely, due to the effective pressure supplied by the 
random velocities that the particles would inevitably ac- 
quire during the collapse when deviations from strict 
spherical symmetry are allowed. From the virial theo- 
rem, energy redistributes among the dust particles such 
that the total kinetic energy becomes —1/2 the total po- 
tential energy in the final state. This suggests that the 
virial radius R v be 1/2 the turn-around radius in the 
matter-dominated universe. Therefore, between turn- 
around and virialization, the density of the sphere in- 
creases by 8. Meanwhile, since the time doubles, the 
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(M 2 p v 



1/3 



(A.2) 



background density decreases by a factor of 4. At turn- 
around, the density contrast is p/p = 9n 2 /16, so the 
virial density becomes p v = 18TT 2 p m . 

Other virial quantities can then be determined in terms 
of the virial density and the total mass of the collapsed 
cloud. Ignoring order-unity numerical factors, the virial 
radius for a mass M is R v = (Af/p^) 1 / 3 , and the virial 
velocity of matter particles is V v = (GM/Ry) 1 ^ 2 = 
(G 3 M 2 p v ) 1 / 6 . Finally, the virial temperature for a given 
mass is 

_ pm T r 2 ^m,G 

v ~ 2k v ~ 2k 
where p m rs 0.6m p is the (appropriately averaged) molec- 
ular mass in the cloud. These are all time dependent 
through the evolution of the background densities, and 
the last relationship gives rise to ([6]). 

A slight generalization of the above argument allows 
the calculation of the virial density in other backgrounds. 
For instance, the same collapsing solution can be used to 
determine the virial density when the background uni- 
verse has curvature; however, the virial density is deter- 
mined by comparison to the spatially curved background. 
An analytic approximation to the exact virial density is 
given in (29|, based on a derivation in [3(3]. In the case 
of a background universe with a cosmological constant 
(but no spatial curvature), the evolution of the overdense 
sphere must be calculated numerically, but the virial den- 
sity is derived in much the same way. One caveat is 
that the cosmological constant modifies the gravitational 
potential, so the virial radius is no longer half of the 
maximum radius of the overdense cloud. Based on cal- 
culations in [3lj |. [2^ | found an approximate analytical 
formula for the virial density in terms of the background 
densities. We have calculated the virial density in terms 
of the background densities for a universe with three com- 
ponents, matter, curvature, and A. 

In all the cases, however, the virial density is within 
two orders of magnitude of the dominant component of 
the background density. Since the virial density only 
enters our calculations through a fractional power, the 
precise value is not terribly important, and we simply 
use the result for the matter-dominated universe. (We 
plan to correct this detail in future work.) 



2. Press-Schechter Collapsed Fraction 

The Press-Schechter collapsed fraction [25| is the frac- 
tion of matter bound in halos of mass M or less, which 
one can approximate as the fraction of density fluctua- 
tions at length scale R with density contrast greater than 
the critical density 6 C . It is given by 

" ' (A.3) 



F(M, t) = erf ( - 

\V2a(M,t), 

where a is the root-mean-square density fluctuation with 
a wavelength given by the comoving radius associated 
with the mass scale M: R 3 = (3/4ir)(M/a 3 p m ). 



Therefore another quantity needed for determining the 
SFR is the linear growth of fluctuations in the argument 
of the Press-Schechter function. In linear theory, the vari- 
ance in overdensity in spheres that enclose a mass M is 



a 2 (M,a) = D{af 



dk 
2^ 



k 2 P(k) 



3h(kR) 



kR 



(A.4) 



This form factorizes the variance into an overall lin- 
earized growth factor and a matter power spectrum, 
which is valid after the end of radiation domination. 
It does not work during radiation domination because 
modes inside and outside the horizon evolve differently 
during that time. After the radiation era, P takes into ac- 
count the different evolutions. There are additional small 
effects after radiation domination but before decoupling 
that we will neglect here. 
One can then integrate 



dk 
2rr 



2 k 2 P{k) 



3ji(fcfl) 
kR 



Q 2 Z(p) 2 



(A.5) 



where the dimensionless parameter p = £ 2 M/Mp is pro- 
portional to the mass in terms of the horizon mass at 
matter-radiation equality fl3| . Here £ = p m /nj is a time- 
independent measure (apart from decays which increase 
the density of photons) of the amount of nonrelativis- 
tic matter in the universe. Empirically, the spectrum is 
found to be well approximated by [l3| 



em 



(9.lp- 2 
(50.5 log (834 + p~ 1/3 ) -92 



-| -1/0.27 

(A.6) 



and Q = dp/ p is the amplitude of fluctuations at hori- 
zon entry, which coincides with the value at the end of 
inflation. 

Next, we compute the linear growth factor, assuming 
that matter dominates over both curvature and A during 
the radiation era. If there were a single overall growth 
factor during the radiation era, we would have for sub- 
horizon scales at sufficiently early times that a < a eq 



D(a) 



1 



3 a 
2a~ 



(A.7) 



eq 



This is a well-known result (see for example [28]). After 
the radiation dominated era (a > a eq ), (|A.7j) generalizes 

to [H 



(A. 



where the spectrum T,(p) accounts for the differential 
growth of different wavelengths during the radiation era. 
Shortly after the end of the radiation era, the first term 
becomes unimportant. The well-known result for matter, 
curvature, and/or A-dominated universes is 



G(t) = H(t) 



a(t) 



da 



(aH) 3 



(A.9) 
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as long as the universe is not pure de Sitter [2g] or dom- 
inated by a form of dark energy that is not a simple 
cosmological constant. The authors of [HI obtain an an- 
alytic approximation for (|A.9[) for their SFR, but that 
approximation is not valid for values of A much different 
that in our universe. Therefore we evaluate the integral 
numerically. 

Nevertheless, it is enlightening to have an analytic for- 
mula to get some intuition for parametric dependences. 
Ref. [l3[ finds that the growth factor for a flat universe 
is approximated by 



D(t) 



M/3 

n 2- x 1/3 

u.z 1/ J, 

Pa 



1 



1.44 



-l/3a 



(A.IO) 



where a = 0.795, pa is the vacuum energy density, and 
x = pa/ p m (t) serves as a dimensionless "time" variable. 

3. Cooling and Star Formation Efficiency 

The other key factor in the HS analysis is the star 
formation rate in a collapsed halo of mass M, denoted 
by s(M,t). HS propose a simple model for gas cooling, 
which matches well with numerical studies of a detailed 
multi-phase model [26| . They find that the star formation 
rate is proportional to the radiative cooling rate. This 
model ignores metal-line cooling and the slow change of 
elemental abundances due to star burning. For our pur- 
poses, metal-line cooling is sufficiently small to neglect, 
while the change of abundances (nuclear fuel depletion) 
is important only for special values of the cosmological 
parameters, so we will ignore both of these effects in this 
paper. This also neglects cooling through H 2 molecules, 
a mechanism which may be important in forming the 
first stars in our universe but is believed to be currently 
negligible. 

From their numerical studies, HS argue for a factoriza- 
tion 

( 0, T„<10 4 K 
s{M,t) = q(t) x I 1/3, 10 4 K < T v < 10 6 5 K (A.ll) 
[ 1, 10 65 K<T V 

and calculate the cooling rate at 10 7 K. Based on numer- 
ical results and physical reasoning about winds in star 
forming regions, HS assume a simple form for the pro- 
portionality factor, 



q(t) cx min hp<?(*) 



(A.12) 



The idea is that at high densities, the star formation 
rate saturates at (slightly less than or equal to) the gas 
consumption time scale t*, approximately 2.1 Gyr. Then 
q is the star formation rate due to radiative processes, 
calculated without regard to this upper limit. 

To compute the cooling rate, consider the cooling time 
at a radius r from the center of the cloud. This is 



t(r) 



-kT 



Pb(r) 



1 



n%n(T) 



(A.13) 



where pb is the baryon density within the cloud, nn is the 
hydrogen number density, and f2(T) is the cooling func- 
tion. We can understand this formula as follows: 3fcT/2 
is the energy per particle, Pb/p is the particle number 
per volume, n 2 H is proportional to the rate of hydrogen 
collisions, and Q(T) is proportional to the energy loss per 
collision. To phrase this in terms of more fundamental 
parameters, we write nn — Xpb/rrih, where X is the hy- 
drogen mass fraction. Following HS, we take the density 
of the cloud to be a power law 



Pb(r) 



(3 - V )f b M 



3-7} 



(A.14) 



where fb = Mb/M is the baryon mass fraction. HS note 
that 77 = 2 is the prediction for a thermal distribution, 
but matching to numerical results gives 1.5 < ?; < 2 with 
77 = 1.65 as the best fit value. 

From the above, we can show that 



dr 



1 



cool 



dt/dr 



1 r 



(A.15) 



Moreover, HS argue that the cooling time should be the 
natural timescale for the gas to virialize, t = R v /V v . 
Combining these results leads to 

(A ' 16) 

For shorthand, we have written 

/(T) 3 kTm% 



(A.17) 



2pX 2 n(T) ' 

From 32], the cooling function is known to be 

f7(10 7 K) = 10" 23 erg-cm 3 ■ s" 1 (A.18) 

with little variation over several orders of magnitude in 
the temperature. 

It is not, however, necessary for our purposes to com- 
pute the cooling and star formation rates to this level 
of detail. The cooling rate (|A.16[) depends on time only 
through the virial density p Vl which is, to within an order 
of magnitude, well-approximated by lOMpH 2 . There- 
fore, it is reasonable to fit the star formation rate (|A.12|) 
by a function of H, which HS give as 



q(t) cx # 3 /" 



fJf2m/3 _|_ Jj2m/3^9/2mi-) 



(A.19) 



with m ~ 6 and H ~ 10-ffo providing good fits to sim- 
ulations. This can be rewritten in the form ([5]). Again, 
we hope to correct this inaccuracy in future work. 



4. Synthesis 

The final expression for the HS star formation rate can 
be simplified by the approximation in terms of elemen- 
tary functions of the error function which appears in the 
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Press-Schechter theory, 



erf(y) »1- 



1 



(A.20) 



1 + ' 

Using constants 5 C = 1.6868 and a = 0.707, 6 HS find that 

(A.21) 

where cr„ stands for fluctuations at the mass scale which 
virializes at temperature 10™ K. The relation between 
mass and temperature at time t is found to be 



M = 



(2fcT//; m ) 3 / 2 ^ (2fc7> m ) 3 / 2 
Pv ~ 10GH[t) 



(A.22) 



where A: is Boltzmann's constant and M depends on the 
Hubble parameter at the given time. Using the measured 
value £/M p = 3.3 x 1CP 28 for the matter abundance, 
the dimensionless mass parameter \x n corresponding to a 

6 The constant a is introduced to modify the Press-Schechter frac- 
tion along the lines of the Sheth-Torman fraction. 



virialization temperature T = 10™ K is found to be 

M „ = 1.1 x 10 3+15 » (Jj^j . (A.23) 



Using (|A.6IA.9[) . the corresponding fluctuation amplitude 
is given by 



a n = Q^ n )D(t) 



(A.24) 



A simpler version of the SFR is also suggested by HS, 
namely 



/j* = q(t) 1 - erf 



a 5 C 
2 



(A.25) 



Our numerical comparisons indicate that this differs very 
little from the version (|A.21[) , so we adopt (|A.25[) for our 
subsequent analysis. This is our final SFR, given earlier 
in 0. 
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